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Abstract —Swarm and modular robotics are an emerging 
area in control of autonomous systems. However, coordinating 
a large group of interacting autonomous agents requires careful 
consideration of the logistical issues involved. In particular, 
inter-agent communication generally involves time delay, and 
bandwidth restrictions limit the number of neighbors with 
which each agent in the swarm can communicate. In this 
paper, we analyze coherent pattern dynamics of groups of 
delay-coupled agents, where the communication network is an 
Erdos-Renyi graph. We show that overall motion patterns for 
a globally-coupled swarm persist under decreasing network 
connectivity, and derive the bifurcation structure scaling re¬ 
lations for the emergence of different swarming behaviors as 
a function of the average network degree. We show excellent 
agreement between the theoretical scaling results and numerical 
simulations. 

The emergence of complex dynamical behaviors from 
simple rules of pairwise interaction in aggregates of agents 
is a remarkable yet widespread phenomenon that appears in 
a multitude of application domains. In biology, aggregates 
form at all spatial and temporal scales, from the microscopic 
(ex., aggregates of bacterial cells or the collective motion 
of skin cells in wound healing) [1], [2], [3] to large-scale 
systems of fish, birds, and even humans [4], [5], [6], In 
control systems, the emerging capabilities of swarms can al¬ 
low groups of relatively small, inexpensive agents to achieve 
tasks that are beyond any individual agent. For example, 
aggregates of locally interacting agents have been proposed 
as a means to create scalable sensor arrays for surveillance 
and exploration [7], [8], [9], [10], [11], [12]; distributed 
sensing [13], [14]; cooperative construction [15], [16]; and 
the formation of reconhgurable modal systems [17], [18], 
[19]. 

The application of multi-agent systems for various tasks 
occurs in parallel with a modeling effort aimed at understand¬ 
ing the emergent properties of swarms. Existing literature 
on the subject provides a great selection of both individual- 
based [5], [6], [20], [4] and continuum models [21], [22], 
[2] of group motion. Many biologically-inspired models are 
based on “zones” of attraction, repulsion, and/or alignment 
interactions between agents [23], [24], [25]. Setting different 
ranges for attraction/repulsion interactions between all agents 
in the swarm or introducing heterogeneous dynamics can 
lead to interesting behaviors, such as splitting of a large 
swarm into smaller groups [26], [27], [28]. Swarming with 
nearest-neighbor alignment is studied numerically in [29]. 
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Stochastic interactions between agents are modeled in [30]. 

Even so, most swarming models do not lend themselves 
to analysis, and there is still limited understanding of how 
group behaviors and structures arise out of simple inter¬ 
agent interactions. However, some progress has been made 
in understanding how the structure and parameters of a 
given swarm model contribute to the aggregate motion of the 
group. For example, [31] presents a simulation-based analy¬ 
sis of the different kinds of motion in a fish-schooling model; 
the authors map phase transitions between different aggregate 
behaviors as a function of group size and maximum number 
of neighbors that influence the motion of each “fish”. In [6], 
the authors use simulations to study transitions in aggregate 
motions of prey in response to a predator attack. 

In this paper, we are motivated by the idea of using a 
swarm of interacting autonomous aerial robots to conduct 
surveillance/monitoring of a specified environment, where 
a human operator uses a high-level control to guide the 
swarm as a whole (e.g., to set the monitoring region), 
while individual agent trajectories are governed by swarm 
interactions. The idea is similar to the reduced-order swarm 
control described in [32], [33] and references therein. Our 
goal is to rigorously characterize the swarm motion patterns 
under a simple but general swarming model, as a function 
of model parameters, so as to exploit them for parametric 
control of the system. We incorporate two key model modifi¬ 
cations for real-world applicability: time delay and restricted 
communication bandwidth. 

Systems of interacting individual agents, whether natural 
and engineered, involve some degree of communication 
delay [34], [35], [36], [37]. Time delay can have significant 
impact on system dynamics, leading to destabilization or syn¬ 
chronization of coupled systems [38], [39], [40]. As shown 
in our earlier work with globally delay-coupled swarms 
of homogeneous [41], [42], [43] and heterogeneous agents 
[44], communication delay can cause emergence of new 
collective motion patterns and, in the presence of noise, lead 
to switching between bistable patterns; this, in turn, can lead 
to instability in robotic swarming systems [25], [45], 

Because global coupling is easier to analyze and a reason¬ 
able assumption in cases of high-bandwidth communication 
and when the number of agents is small, many models 
make the mathematically simple but physically implausible 
assumption that swarms are globally coupled (that is, each 
agent is influenced by the motion of all other agents in the 
swarm) [46], [47], [6], [48], However, for large groups of 
agents, bandwidth restrictions generally mean that it is not 
feasible to maintain all-to-all communication. In this pa- 



per, we generalize our previous results for globally-coupled 
swarms by assuming an Erdos-Renyi communication net¬ 
work, where the mean number of neighbors for each agent 
in the swarm is fixed. 

We use mean-field dynamics to analytically predict transi¬ 
tions between regimes of different collective swarm motions 
as a function of model parameters for swarms consisting of 
homogeneous delay-coupled agents with a fixed undirected 
Erdos-Renyi communication network. We use a simple yet 
general particle swarming model, which combines agent 
self-propulsion and inter-agent attraction. This choice can 
be justified in the case of swarming robots with very fast 
relaxation times, such as quadcopters, which can be treated 
as holonomic vehicles over the spatio-temporal scales re¬ 
quired for a surveillance/monitoring deployment. We show 
that the stable motion patterns observed in the globally- 
coupled system (translation, ring formation, and rotation 
about a common center of mass) [42] persist under non- 
global coupling, for appropriately chosen model parameters. 
Our results are verified through numerical simulations. 

I. Problem Statement 

Consider a swarm of delay-coupled agents in the plane, 
with each agent indexed by i. We build an Erdos-Renyi 
communication network on the swarm as follows. Starting 
with a globally-coupled network, pick an existing link at 
random, using a uniform distribution over the existing links, 
and remove it, until the mean degree matches some target 
value. Let AT t denote the set of neighbors of agent i, that is, 
the set of agents with which agent i shares a communication 
link. Because of communication/sensing delays, information 
about the position of the neighboring agents is available 
with delay r, assumed to be the same for all agents. For 
simplicity, we assume that all agents in the swarm have 
identical dynamics and that all connections in the network 
are bidirectional and fixed in time. The agents have self¬ 
propulsion and are attracted to their neighbors with strength 
that depends on the coupling coefficient a. Then the motion 
of agent i is governed by the following equation: 

*» = (!- N 2 )i-i - X ( r *(*) “ r f W)> (!) 

feW 

where superscript r is used to denote time delay, so that 
rj (t) = Tj(t — t), and | • | denotes the Euclidean norm. 
The first term in the above equations represents the self¬ 
propulsion of swarm agents. The second term models pair¬ 
wise attraction between each agent and its neighbors in the 
swarm. This simplified model does not include short-range 
repulsion or other collision-avoidance strategies; however, 
earlier studies indicate that the collective dynamics of the 
swarm are not significantly altered by the introduction of 
short-range repulsion terms [42], 

We examine the dynamics of the system analytically in 
the limit where the system is almost completely connected 
( ( Ar ~ 1 )~™ rd W») x, where card(Ni) is the number of 
neighbors of node i ), and show via simulations that the 
approximations made in the almost-connected limit hold 


closely even as the mean coupling degree is reduced to less 
than 50% of possible links. 

II. System Dynamics in the Mean-Field 

In [42] we derived the mean-field dynamics of the system 
in the limit N —» oo where the delayed coupling was 
considered to be all-to-all. The difference considered in the 
model given by Eq. [I] is that each agent’s motion depends 
only on its neighbors, rather than all other agents in the 
swarm. Let R(t) = Y^iLi r *(^) d en °te the position of the 
center of mass of the swarm, and let Sr.i(t) = rj(f) — R(f); 
note that = 0- Let (■, •) denote the dot product in 

R 2 . The motion of the center of mass is governed by 

R = (1- |R| 2 )R 

1 N 

- n E + (l ^ 2 + 2(R ’ + "*)) 

V i= 1 

- f ( R - Rr) - ^ E - E ) * 

i=i V feW ) 
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where pi = card(Mi ) is the number of neighbors of 
agent i (in the globally coupled case, pi = N — 1) and 
p = Pi ' s the mean degree of the network. We 

denote the mean fraction of missing links by e, so that 
p = (1 — e)(N — 1). Note that, if pt = p is the same for all 
i, 

N / \ N N 

X !\ p6ri ~ X d ' r J = - E E St j = -pX 6r i = 0 
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since Xa=i hif: = 0 and in the double summation, each hr, 
appears p times because the network is undirected and each 
node has p neighbors. However, in our network p, is drawn 
from a binomial distribution. Our previous work shows that 
in many instances either individual deviations from the center 
of mass are small, or in aggregate they average out over 
the whole population. (We discuss situations in which this 
assumption breaks in a later section). Then, neglecting all 
terms of order hr, and in the limit N —> oo, the center of 
mass motion is given approximately by 

R = (1 - |R|)R - a( 1 - e)(R(t) - R T (i)). (4) 

We note that the mean-field dynamics are identical to those 
for the globally-coupled homogeneous swarm described in 
[42], with coupling constant a replaced by an ‘effective’ 
value a e ff = al l - e). Quite remarkably, simulation results 
indicate that the system exhibits similar collective motions 
to the globally coupled case, even as p is significantly 
decreased. These basic collective motion patterns are “trans¬ 
lation”, where the entire swarm collapses to a point and 
travels along a straight-line trajectory at constant speed; 
“ring” motion, where the swarm agents form concentric 
counter-rotating rings about the stationary center of mass; 
and “rotating” motion, where the agents collapse to a small 




volume and collectively rotate about a fixed point. The 
collective motions of the swarm and the effects of non- 
global coupling are described in more detail in the following 
section. 


III. Collective swarm motions 

As in the globally coupled case, the steady-state motions 
of the swarm depend on values of the coupling coefficient a, 
the delay r; in addition, there is now a dependence on the 
fraction of missing links e. The collective motion patterns of 
the swarm for different values of the parameters a and r are 
described in more detail below. 


A. Translating state 

In the translating state, the agent locations all lie close to 
the swarm center of mass, and the swarm moves with con¬ 
stant speed and direction. Following the calculation in [42], 
it can be shown that the translation speed ||R|| must satisfy 
|| R|| 2 = 1— a e ffT. The system exhibits a pitchfork bifurcation 
at a e ffT = 1, where the translating state disappears. 


B. Ring state 

For all values of a and r, 0 admits a stationary solution, 
R(i) = R(0). In this state, the agents converge to a pair of 
concentric, counter-rotating rings about the stationary center 
of mass. The mean radius and angular velocity of the agents 
in the ring state satisfy 


P=l/\/^i (5) 

w = ±\/orff, (6) 


respectively (see Fig. [2]). The stability of the ring state is 
determined by the characteristic equation 


M( A; a, r) 


(A 2 —A+a e ff ) 2 —(A 2 —A+a e ff)a e ffe Ar 


(7) 


and is lost along Hopf bifurcation curves which give rise to 
a rotating state. 


C. Rotating state 

Solving for values where roots of M cross the imaginary 
axis, we obtain Hopf bifurcation curves in the a and r plane: 


1 

v / 2a e ff — 1 


^arctan 


x/2a e ff — 1 
1 — «eff 


2m7r 


( 8 ) 


to £ Z. When e = 0, we recover the equations for the 
globally coupled system by taking a e ff = a; the factor of 
(1 — e) in a e ff which results from breaking a fraction of the 
links in the global network represents a perturbation from the 
globally-coupled case (Fig. [lj. Note that the pitchfork and 
Hopf bifurcation curves meet at a Bogdanov-Takens point 
when a = r = 2(1 — e). 

The rotating state, in which the agents collapse and 
collectively rotate about a fixed point is created at the first 
Hopf bifurcation. In the case of global coupling, all agent 
positions coincide; however, when coupling is not global, 
different agents circle the fixed point with equal angular 
frequency but have different radii, and have a fixed relative 
phase offset from the center of mass, depending on their 



Fig. 1: Hopf bifurcation curves for e = 0 (blue), e = 0.1 
(green), and e = 0.2 (red). The dashed lines show the 
location of the pitchfork bifurcation where the translating 
state disappears. In the mean field, the translating state occurs 
in region I; the ring state in region II, and the rotating state, 
in region III. 


coupling degree (see Fig. [3j. The radius and angular velocity 
of the center of mass of the swarm in the rotating state satisfy 

1 f sinwr 

P = TA V 1 ~ °e£f- (9) 

\uj\ V lj 

w 2 = a e ff(l - coswr) (10) 

(see Fig. |4|. 

We now investigate the precise spatial organization of 
swarming agents in the rotating state. Using rj = R r + <5rJ 
allows us to write the coupling of agent i to the other 
particles as 

~ RT ) - SrT j 

j&Mi jeMi 

~ Pi[ r i — R T ), (11) 

after neglecting the sum of small terms <5rJ. In this way, the 
dynamics of agent i becomes 

ri = (1 - N 2 )U - - RT W), (12) 

so that dynamically, less-than-global coupling in our Erdos- 
Renyi network is approximately equivalent to a fully con¬ 
nected network of agents with heterogeneous coupling coef¬ 
ficients 

ai = api/{N- 1), (13) 

where the factor of N - 1 is introduced so that a,i = a in 
the fully connected case where pi = N — 1. In the limit 
N —> oo, we have 

Vi = (1 - |fj| 2 )rj - at (r, - R r ), (14) 

where the motion of the center of mass R is given by (|4j. 

Let (p,9) and ( pi,9i ) denote the polar coordinates of the 
swarm center of mass and of agent i, respectively. In these 
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Fig. 2: Comparison of the theoretical values of radius (top) 
and angular velocity (bottom) of the agents in the ring 
state, for different connection degrees, compared with full- 
swarm simulation of 150 agents. The x-axis represents the 
percentage of active links, out of all possible links; all links 
are bidirectional and time-invariant. Error bars are shown 
one standard deviation above and below the mean values. 
The simulations were run with a = 1 and r = 2.5. 


coordinates, the equations of motion for the swarm are 


p = (1 — p 2 0 2 — p 2 )p + p8 2 — a(p — p T cos (9 — 6 T )) 

(15) 

p9 = (1 — p 2 0 2 — p 2 )p9 — 2 p9 — ap T sin(0 — 9 T ) (16) 

Pi = (1 - Pi^i ~ Pi)Pi + Prfi - a i(Pi - P T cos (0i - 6 T )) 

(17) 

PiOi = (1 - p 2 9 2 - p 2 )pidi - 2 piOi - dip T sin(0i - 9 T ), 

(18) 


where a = YliLi a i I s th e mean coupling coefficient. 

Without loss of generality, we set the origin of the polar 
coordinate system at the center of rotation so that p and pi 
are constant in time. Also, we have 9 = 0j = 0. Let ui = 9 
and uii = 8i. For convenience, we define £* = a,/ a and 
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Fig. 3: Snapshot of simulation showing swarm in rotating 
state about a stationary center point (marked by ‘x’). The 
colors indicate the coupling degree of each agent, normalized 
by iV — 1. Here a = 1, r = 4.5, and p/(N — 1) = 0.75. 


C P z Pi / P- The equations of motion then become 

0 = pu 2 — ap(l — cos(wr)) (19) 

0 = (1 — p 2 uj 2 )pu — apsin(wr) (20) 

0 = CpPul - £a(? p p 

— pcos((oji — uj)t + 9i{ 0) — 0(0) + ujt)) (21) 

o = (i 

— ^dpsin((wj — uj)t + 0j(O) — 0(0) + ujt)). (22) 

Note that equations ( [2T] > and ( |22) > can only be satisfied for all 
times t if u>i = ui. Let A0^ = 0^(0) —0(0) denote the angular 
offset between particle i and the center of mass. Simplifying, 
we finally have 

0 = w 2 — a(l — cos(wr)) (23) 

0 = (1 — pi 2 uj 2 )uj — asin(wr) (24) 

0 = C p w 2 - C a a(C p - cos(wt + A9i)) (25) 

0 = (1 - (Q 2 p 2 uj 2 )CpUJ - C a asin(ujT + A0,)). (26) 

This set of coupled nonlinear equations can be solved numer¬ 
ically for p, u>, , and A0, for different values of d, r and 

£®. We compare numerical solutions for the swarm center 
of mass radius and angular velocity with simulations of a 
swarm with Erdos-Renyi communication network structure 
(see Fig. |dji. Solution curves for £* and A0, for different 
values of a are shown in Fig. [5j a comparison with sim¬ 
ulation results is shown in Fig7[6] The slight discrepancy 
in the rotating state radius in Fig. |d] and in £* in Fig. |6] 
is understood as follows. The radius of the center of mass 
computed from ( |23] > and ( |24| ) assumes that agent positions’ 
deviate only slightly from the center of mass. However, as 
the mean coupling coefficient decreases the agents become 
spread out over an extended arc (as seen in Fig. [3]) and 
our assumption becomes invalid. In this ‘arc’ configuration 
the center of mass of the swarm moves closer towards the 
center of rotation than theory predicts. The analogue to a 
system with perturbed coupling coefficient breaks down here; 





















for a globally-connected swarm with decreasing coupling 
coefficient a, the rotating state disappears when the system 
crosses the curve dr 2 = 2, where the rotating state radius 
diverges. The swarm then transitions to a translating state. It 
is, however, remarkable, that the mean-field analysis captures 
so much of the overall swarm behavior even as the coupling 
degree is significantly decreased. 
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Fig. 4: Comparison of the theoretical values of radius (top) 
and angular velocity (bottom) of the agents in the rotating 
state with full-swarm simulations of 150 agents, for different 
connection degrees. The simulation values are obtained by 
averaging over all agents int he swarm. The x-axis represents 
the percentage of active links, out of all possible links; all 
links are bidirectional. Error bars are shown one standard 
deviation above and below the mean. 


IV. Conclusion 

In this paper we have analyzed the collective motion 
patterns of a swarm with Erdos-Renyi communication net¬ 
work structure, using a mean-field approach from statistical 
physics, with the assumption that the number of agents goes 
to infinity. We derived bifurcation diagrams demarcating 
regions of different collective motions, for different values of 
mean degree in the communication network. We showed that 
behaviors described in [42] for the globally-coupled swarm, 
namely translation, ring state, and rotation, persist as commu¬ 
nication links are broken, even though the bifurcation curves 
are shifted as coupling degree of the network decreases. 




Fig. 5 : Theoretical values for ratio of radius for agent i to 
radius of swarm center of mass (top) and phase difference 
(bottom) as a function of £* = dj/d, for a = 1, 0.8, 0.6, 
and 0.4. 


We derive expressions for the speed of the swarm in the 
translating state as a function of time delay and coupling 
coefficient; for the mean radius and angular velocity of 
agents in the ring state; and for the angular velocity, and 
individual radii and phase offsets for individual agents in 
the rotating state. We have verified these calculations with 
simulations of the full-swarm dynamics. It is remarkable 
that our model reduction, which starts with N second-order 
delay-differential equations and yields one equation of the 
same type, is able to quantitatively capture so many aspects 
of the full swarm dynamics, even as the coupling degree of 
agents within the swarm is significantly decreased. 

In the case that many agents are coordinating together, 
limited communication bandwidth makes all-to-all commu¬ 
nication infeasible, and may lead to significant communi¬ 
cation delays. By dropping the requirement for all-to-all 
communication used in our previous work, the current paper 
brings us one step closer to a possible implementation of 
swarming control algorithms for very large aggregates of 
agents. Understanding the natural emerging dynamics of the 
system in these circumstances allows us to exploit them when 
designing controls for swarming applications. 

In future work, we will study the effects of using nearest- 
neighbor communication rather than the random Erdos-Renyi 
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Fig. 6: Comparison of simulation and theoretical values for 
radius (top) and phase difference from the center of mass 
(bottom) for agent i. Results are shown for two swarms, 
the first with a = 0.6 (left) and the second with a = 0.4 
(right) (that is, 60% of links active and 40% of links active, 
respectively). For both swarms, r = 4.5 sec, and number of 
agents is 150. All links are bidirectional. Note that within 
each swarm, agents with higher coupling degree lie further 
from the center of rotation. 


communication network used here. We will verify our results 
in the lab using a swarm of quadcopters operating in a mixed 
real-virtual environment, in which arbitrarily large numbers 
of virtual agents can be simulated to interact with the group 
of real vehicles. We will test how our results scale with 
the number of agents in the network, and apply parametric 
control for dynamic pattern-switching. 

Acknowledgments 

This research was performed while KS held a National 
Research Council Research Associateship Award at the U.S. 
Naval Research Laboratory. This research is funded by the 
Office of Naval Research contract no. N0001412WX2003 
and the Naval Research Laboratory 6.1 program contract no. 
N0001412WX30002. LMR is a post-doctoral fellow at Johns 
Hopkins University supported by the National Institutes of 
Health. 

References 

[1] E. O. Budrene and H. C. Berg, “Dynamics of formation of symmetrical 
patterns by chemotactic bacteria,” Nature, vol. 376, pp. 49-53, 1995. 

[2] A. A. Polezhaev, R. A. Pashkov, A. I. Lobanov, and I. B. Petrov, 
“Spatial patterns formed by chemotactic bacteria Escherichia coli,” The 
International Journal of Developmental Biology, vol. 50, pp. 309-314, 
Jan. 2006. 

[3] R. M. Lee, D. H. Kelley, K. N. Nordstrom, N. T. Ouellette, and 
W. Losert, “Quantifying stretching and rearrangement in epithelial 
sheet migration.” New Journal of Physics, vol. 15, Feb. 2013. 

[4] K. r. Tunstrp m, Y. Katz, C. C. Ioannou, C. Huepe, M. J. Lutz, and 
I. D. Couzin, “Collective states, multistability and transitional behavior 
in schooling fish,” PLoS computational biology, vol. 9, p. el002915, 
Jan. 2013. 

[5] D. Helbing and P. Molnar, “Social force model for pedestrian dynam¬ 
ics,” Physical Review E, vol. 51, no. 5, pp. 4282^1286, 1995. 


[6] S.-H. Lee, “Predator’s attack-induced phase-like transition in prey 
flock,” Physics Letters A, vol. 357, pp. 270-274, Sept. 2006. 

[7] P. Bhatta, E. Fiorelli, F. Lekien, N. E. Leonard, D. A. Paley, F. Zhang, 
R. Bachmayer, D. M. Fratantoni, R. E. Davis, and R. J. Sepulchre, 
“Coordination of an underwater glider fleet for adaptive sampling,” in 
Proceedings of the International Workshop on Underwater Robotics, 
no. August, pp. 61-69, 2005. 

[8] Y.-l. Chuang, M. R. DOrsogna, D. Marthaler, A. L. Bertozzi, and L. S. 
Chayes, “State transitions and the continuum limit for a 2D interacting, 
self-propelled particle system,” Physica D: Nonlinear Phenomena, 
vol. 232, pp. 33^17, Aug. 2007. 

[9] K. M. Lynch, I. B. Schwartz, P. Yang, and R. A. Freeman, “Decen¬ 
tralized Environmental Modeling by Mobile Sensor Networks,” IEEE 
Transactions on Robotics, vol. 24, pp. 710-724, June 2008. 

[10] P. Yang, R. A. Freeman, and K. M. Lynch, “Multi-agent coordina¬ 
tion by decentralized estimation and control,” IEEE Transactions on 
Automatic Control, vol. 53, no. 11, pp. 2480-2496, 2008. 

[11] W. Wu and F. Zhang, “Cooperative exploration of level surfaces of 
three dimensional scalar fields,” Automatica, vol. 47, pp. 2044-2051, 
Sept. 2011. 

[12] W. Wu and F. Zhang, “Robust Cooperative Exploration With a 
Switching Strategy,” IEEE Transactions on Robotics, vol. 28, no. 4, 
pp. 828-839, 2012. 

[13] T. H. Chung, J. W. Burdick, and R. M. Murray, “A decentralized 
motion coordination strategy for dynamic target tracking,” in IEEE 
International Conference on Robotics and Automation, vol. 2006, 
pp. 2416-2422, 2006. 

[14] S. Kar and J. M. F. Moura, “Distributed linear parameter estimation 
in sensor networks: Convergence properties,” in 2008 42nd Asilomar 
Conference on Signals, Systems and Computers, pp. 1347-1351, Ieee, 
Oct. 2008. 

[15] F. Augugliaro, S. Lupashin, M. Hamer, C. Male, M. Hehn, M. W. 
Mueller, J. S. Willmann, F. Gramazio, M. Kohler, and R. D’Andrea, 
“The flight assembled architecture installation: Cooperative conduc¬ 
tion with flying machines,” IEEE Control Systems, vol. 34, no. August, 
pp. 46-64, 2014. 

[16] J. Werfel, K. Petersen, and R. Nagpal, “Designing collective behavior 
in a termite-inspired robot construction team,” Science (New York, 
N.Y.), vol. 343, pp. 754-8, Feb. 2014. 

[17] S. Bandyopadhyay, S.-j. Chung, and F. Y. Hadaegh, “Probabilistic 
Swarm Guidance using Inhomogeneous Markov Chains,” 2014. 

[18] M. Dorigo, D. Floreano, L. M. Gambardella, F. Mondada, S. Nolfi, 
T. Baaboura, M. Birattari, M. Bonani, M. Brambilla, A. Brutschy, 
D. Burnier, A. Campo, A. L. Christensen, A. Decugniere, G. Di 
Caro, F. Ducatelle, E. Ferrante, A. Forster, J. M. Gonzales, J. Guzzi, 
V. Longchamp, S. Magnenat, N. Mathews, M. Montes de Oca, 
R. O’Grady, C. Pinciroli, G. Pini, P. Retornaz, J. Roberts, V. Sperati, 
T. Stirling, A. Stranieri, T. Stutzle, V. Trianni, E. Tuci, A. E. Turgut, 
and F. Vaussard, “Swarmanoid: A Novel Concept for the Study 
of Heterogeneous Robotic Swarms,” IEEE Robotics & Automation 
Magazine, vol. 20, pp. 60-71, Dec. 2013. 

[19] M. Rubenstein, A. Cornejo, and R. Nagpal, “Programmable self- 
assembly in a thousand-robot swarm,” Science, vol. 345, pp. 795-799, 
Aug. 2014. 

[20] T. Vicsek, A. Czirok, E. Ben-Jacob, I. Cohen, and O. Shochet, “Novel 
type of phase transition in a system of self-driven particles,” 2006. 

[21] L. Edelstein-Keshet, D. Grunbaum, and J. Watmough, “Do travelling 
band solutions describe cohesive swarms? An investigation for migra¬ 
tory locusts,” Journal of Mathematical Biology, vol. 36, pp. 515-549, 
July 1998. 

[22] C. M. Topaz and A. L. Bertozzi, “Swarming Patterns in a Two- 
Dimensional Kinematic Model for Biological Groups,” SIAM Journal 
on Applied Mathematics, vol. 65, pp. 152-174, Jan. 2004. 

[23] J. M. Miller, A. Kolpas, J. P. J. Neto, and L. F. Rossi, “A continuum 
three-zone model for swarms,” Bulletin of mathematical biology, 
vol. 74, pp. 536-61, Mar. 2012. 

[24] I. Tarras, N. Moussa, M. Mazroui, Y. Boughaleb, and A. Hajjaji, 
“Collective Behavior of Interacting Particles: Radius-Dependent Phase 
Transition,” Modern Physics Letters B, vol. 27, p. 1350028, Feb. 2013. 

[25] C. Viragh, G. Vasarhelyi, N. Tarcai, T. Szorenyi, G. Somorjai, 
T. Nepusz, and T. Vicsek, “Flocking algorithm for autonomous flying 
robots,” Bioinspiration & biomimetics, vol. 9, p. 025012, June 2014. 

[26] M. Kumar, D. P. Garg, and V. Kumar, “Segregation of Heterogeneous 
Units in a Swarm of Robotic Agents,” IEEE Transactions on Automatic 
Control, vol. 55, no. 3, pp. 743-748, 2010. 























[27] Z. Chen, T. Chu, and J. Zhang, “Swarm splitting and multiple targets 
seeking in multi-agent dynamic systems,” 49th IEEE Conference on 
Decision and Control (CDC), pp. 4577-4582, Dec. 2010. 

[28] Z. Chen, H. Liao, and T. Chu, “Clustering in multi-agent swarms 
via medium-range interaction,” EPL (Europhysics Letters), vol. 96, 
p. 40015, Nov. 2011. 

[29] T. Vicsek, A. Czirok, E. Ben-Jacob, I. Cohen, and O. Shochet, 
“Novel Type of Phase Transition in a System of Self-Driven Particles,” 
Physical Review Letters, vol. 75, no. 6, pp. 1226-1229, 1995. 

[30] C. Nilsen, J. Paige, O. Warner, B. Mayhew, R. Sutley, M. Lam, 
A. J. Bemoff, and C. M. Topaz, “Social aggregation in pea aphids: 
experiment and random walk modeling,” PloS ONE, vol. 8, p. e83343, 
Jan. 2013. 

[31] S. V. Viscido, J. K. Parrish, and D. Griinbaum, “The effect of 
population size and number of influential neighbors on the emergent 
properties of fish schools,” Ecological Modelling, vol. 183, pp. 347- 
363, Apr. 2005. 

[32] C. Belta and V. Kumar, “Abstraction and control for groups of robots,” 
IEEE Transactions on Robotics, vol. 20, no. 5, pp. 865-875, 2004. 

[33] N. Michael, C. Belta, and V. Kumar, “Controlling three dimensional 
swarms of robots,” Proceedings - IEEE International Conference on 
Robotics and Automation, vol. 2006, no. May, pp. 964-969, 2006. 

[34] A. Martin and S. Ruan, “Predator-prey models with delay and prey 
harvesting,” Journal of Mathematical Biology, vol. 43, pp. 247-267, 
Sept. 2001. 

[35] S. Bernard, J. Belair, and M. C. Mackey, “Bifurcations in a white- 
blood-cell production model,” Comptes Rendus Biologies, vol. 327, 
pp. 201-210, Mar. 2004. 

[36] N. A. M. Monk, “Oscillatory Expression of Hesl, p53, and NF-kB 
Driven by Transcriptional Time Delays,” Current Biology, vol. 13, 
pp. 1409-1413, 2003. 

[37] E. Forgoston and I. B. Schwartz, “Delay-induced instabilities in self- 
propelling swarms,” Physical Review E, vol. 77, no. 302, p. 035203, 
2008. 

[38] A. Papachristodoulou and A. Jadbabaie, “Synchronization in Oscillator 
Networks with Heterogeneous Delays, Switching Topologies and 
Nonlinear Dynamics,” in Proceedings of the 45th IEEE Conference 
on Decision and Control, pp. 4307^-312, Ieee, 2006. 

[39] Z. Zuo, C. Yang, and Y. Wang, “A unified framework of exponential 
synchronization for complex networks with time-varying delays,” 
Physics Letters A, vol. 374, no. 19-20, pp. 1989-1999, 2010. 

[40] A. Englert, S. Heiligenthal, W. Kinzel, and I. Kanter, “Synchronization 
of chaotic networks with time-delayed couplings: An analytic study,” 
Physical Review E - Statistical, Nonlinear, and Soft Matter Physics, 
vol. 83, no. 4, pp. 1-14, 2011. 

[41] L. Mier-y-Teran Romero, E. Forgoston, and I. B. Schwartz, “Noise, 
Bifurcations, and Modeling of Interacting Particle Systems,” in Pro¬ 
ceedings of the IEEE/RSJ International Conference on Intelligent 
Robots and Systems, pp. 3905-3910, Jan. 2011. 

[42] L. Mier-y-Teran Romero, E. Forgoston, and I. B. Schwartz, “Coher¬ 
ent Pattern Prediction in Swarms of Delay-Coupled Agents,” IEEE 
Transactions on Robotics, vol. 28, pp. 1034-1044, Oct. 2012. 

[43] B. S. Lindley, L. Mier-y-Teran Romero, and I. B. Schwartz, “Noise 
Induced Pattern Switching in Randomly Distributed Delayed Swarms,” 
in Proc. 2013 American Control Conference (ACC 2013), pp. 4587- 
4591, 2013. 

[44] K. Szwaykowska, L. Mier-y-Teran Romero, and I. B. Schwartz, 
“Collective motions of heterogeneous swarms,” IEEE Transactions on 
Automation Science and Engineering (to appear). Sept. 2015. 

[45] Y. Liu, K. M. Passino, and M. Polycarpou, “Stability analysis of one¬ 
dimensional asynchronous swarms,” IEEE Transactions on Automatic 
Control, vol. 48, no. 10, pp. 1848-1854, 2003. 

[46] S. Motsch and E. Tadmor, “A New Model for Self-organized Dynam¬ 
ics and Its Flocking Behavior,” Journal of Statistical Physics, vol. 144, 
pp. 923-947, Aug. 2011. 

[47] Y. Chen and T. Kolokolnikov, “A minimal model of predator-swarm 
interactions.,” Journal of the Royal Society, Interface/the Royal Soci¬ 
ety, vol. 11, p. 20131208, May 2014. 

[48] F. Vecil, P. Lafitte, and J. Rosado Linares, “A numerical study of 
attraction/repulsion collective behavior models: 3D particle analyses 
and ID kinetic simulations,” Physica D, vol. 260, pp. 127-144, 2013. 



